Anisotropic plasmons in a two-dimensional electron gas with spin-orbit interaction 
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Spin-orbit coupling induced anisotropies of plasmon dynamics are investigated in two-dimensional 
semiconductor structures. The interplay of the linear Bychkov-Rashba and Dresselhaus spin-orbit 
interactions drastically affects the plasmon spectrum: the dynamical structure factor exhibits vari- 
ations over several decades, prohibiting plasmon propagation in specific directions. While this 
plasmon filtering makes the presence of spin-orbit coupling in plasmon dynamics observable, it also 
offers a control tool for plasmonic devices. Remarkably, if the strengths of the two interactions are 
equal, not only the anisotropy, but all the traces of the linear spin-orbit coupling in the collective 
response disappear. 
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Spin-orbit coupling in semiconductor heterostructures 
has received wide attention recently - it has been inves- 
tigated as a source of new fundamental spin physics as 
well as a control interaction in spintronics applicationsi'-S. 
Two spin-orbit terms are relevant in zinc-blende systems 
exemplified by two-dimensional GaAs or InAs electron 
gases: the Bychkov-Rashba^ interaction (coupUng con- 
stant a), which is due to the structure inversion asym- 
metry, and the Dresselhaus interaction- (couphng con- 
stant j3), which is due to the bulk inversion asymme- 
try^. Alone, these interactions lead to an isotropic single- 
particle and plasmon spectrum. Taken together, they 
imprint the underlying heterostructure anisotropy onto 
the single and many-particle properties. Most studies of 
the spin-orbit coupling effects have been on the single- 
particle level. While the presence of spin-orbit coupling 
leads to such notorious effects as spin relaxation^ or 
spin Hall currents^"^, fascinating phenomena originate 
from the interplay of the Bychkov-Rashba and Dressel- 
haus terms. The interplay often leads to pronounced 
anisotropies2i2ii2iiiS'iiii^'i^, but this is not a rulc^"'. 

Recently several many-body effects important for spin 
properties of semiconductor nanosystems have been stud- 
ied in 2DESi^ii^iii. One of the key phenomena due 
to spin-orbit interaction (SOI) in many-spin systems is 
the generation of the inter-chirality-subband electron- 
hole continuum. However, the dispersive and dissipative 
modifications, induced by individual (Bychkov-Rashba or 
Dresselhaus) SOI, are difficult to observe in experiment 
- their effect is isotropic and proportional to the small 
SOI couplingi^>iai20^^^.24^_ In real samples the in- 
terplay of different SOI mechanisms takes place and as we 
show here, it results in the striking anisotropy effect on 
the spectral properties of collective excitations in 2DES. 
This qualitatively strong effect can serve as a valuable tool 
to facilitate the observation and exploitation of usually 
weak SOI effects on many-body properties of 2DES. 

An important outcome of our theory is the prediction 
of plasmon directional filtering: the interplay of the spin- 
orbit couplings leads to plasmon overdamping (blocking) 



in certain special directions of propagation and for cer- 
tain magnitudes of the wave vector. This may be surpris- 
ing at first sight, given that the spin-orbit effects on the 
plasmon dispersion and on the electron-hole excitation 
energies are in themselves quite small. However even 
small energy shifts are sufficient, at these special wave 
vectors, to move the plasmon in or out of resonance with 
electron-hole excitations, thus producing a large effect 
on the plasmon damping. By scanning for plasmons in 
different directions, this distinct absence of propagation 
in certain directions should be experimentally verifiable, 
since the dynamical structure factor varies by orders of 
magnitude as a function of the propagation angle. In 
addition to making the spin-orbit presence experimen- 
tally visible, the anisotropy is attractive for plasmonics 
designs as a substitute for surface patterning to achieve 
directional plasmon propagations^. This prospect is en- 
forced by the possibility to control - even turn on and off 
- plasmon propagation: both a and (3 can be tuned by 
external gatesi (see also^2i22,) allowing for the anisotropy 
to be tailored. In fact, the anisotropy vanishes (filtering 
turned off) for a = ±/3. More surprising, in this case 
the (linear) spin-orbit couplings play no role in plasmon 
dynamics - the isotropic contributions by the individual 
spin-orbit terms cancel each other. 

We calculate the effect of joint Bychkov-Rashba and 
Dresselhaus SOI on the propagation of plasmons in the 
(001) plane of a zincblende semiconductor heterostruc- 
ture. We consider samples at low temperatures with high 
density 2DES where the kinetic energy of electrons dom- 
inates the Coulomb potential energy. In this regime it 
is legitimate to neglect the effect of exchange and corre- 
lations in treating plasmon excitations. We use the ran- 
dom phase approximation^^ and calculate the anisotropic 
Lindhard polarization function for a given wave vector 
q and frequency to. The space in which the imaginary 
part of the Lindhard function differs from zero is known 
as the electron-hole continuum (EHC)^*^, for it describes 
the spectrum of electron-hole excitations. The interplay 
of the Bychkov-Rashba and Dresselhaus SOI leads to the 
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appearance of several sub-regions of the EHC separated 
by boundaries across which the imaginary part of the 
dielectric function exhibits sharp variations. An inter- 
esting effect arises when the frequency of a plasmon of a 
given q but variable propagation direction crosses these 
boundaries: The sudden rise in the density of electron- 
hole excitations causes strong Landau damping, actually 
overdamping the plasmons over a range of wave vector 
orientations. This anisotropy of the plasmon spectrum 
should be observable through the pronounced anisotropy 
of the dynamical structure factor, as shown below. 
Our spin-orbit interaction Hamiltonian is^ 

Hsoi = a {cFxky ~ CFykx) + (3 [crxkx - ^yky) , (1) 

where ax,y are the Pauli matrices, k is the in-plane elec- 
tron momentum with magnitude k and polar angle 0k- 
The eigenvectors of the Hamiltonian H = Hq + i?soi 
with Hq = k'^ /2m* (m* is the electron effective mass 
and ^ = 1) are 



(2) 
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corresponding to the single-particle spin-split branches of 
the electron energy, 
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2m* 



(3) 



labeled by the chirality /i = ±1; A is the area of the 
2DEG. The phase of the spinor in Eq. is ^}{a, (3, 0k) = 
Arg[ae*'^'' -I- i/?e~*"^''] and the angle dependent Rashba- 
Dresselhaus momentum is 



e, 0k) = /9v/l + sin(20)sin(20k), 
with amplitude p — m* \J o? 



(4) 



. The angle parameter 
defined as tan^ ~ (3 /a, describes the relative strength 
of the Bychkov-Rashba and Dresselhaus SOI. The Fermi 
momenta of the subbands ^ are also angle dependent: 



fc^^(p, 0, 0k) - V2mi?f +C(p,^,0k)2 - /i 0, <j>^) , 

(5) 

where the total carrier density n determines the Fermi 



energy. 



Figure [T] illustrates 



the energy spectrum of the chirality subbands and the 
anisotropy of the Fermi contour (note that the Fermi en- 
ergy can be negative). 

The Lindhard polarization function^^ in the presence 
of SOI is defined as a sum over chirality indices n(g, uj) — 

J2t,,u=±i'^i^'^i?^^)^ with 



dk f[E^{k)]-f[E,{k + cf)] 



{2nY E^{k) 
^J^^iiy (k,k + q 



E^{k + q) + ijj + i{) 



(6) 



where f[E^{k)] is the Fermi distribution function. The 



form factors 



,y {k, k + qj 



^i^v (k, k + qj 



q ) are given by 
1 



[1 + fj,iy cos (A(/7q)] , 



(7) 
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FIG. 1: Fermi contours in the momentum plane {kx,ky) for a 
spin-orbit interaction of the form given in Eq. ([l]): (a,) a ^ /3 
and (b) a = f3. The arrows indicate the spin direction. 



where we define Atpq = ip{a, /3, 0k) — fia, P, 0k-i-q)- No- 
tice that in contrast to the case of pure Bychkov-Rashba 
or pure Dresselhaus SOI, here the polarization function 
depends both on the magnitude, q, and orientation, 0q, 
of the wave vector q. Making the replacement k —f —k—q 
in the term of ^ with f[E^{k + q)] and regrouping, we 
can represent the polarization function in the compact 
form Il{q,uj) = J2u,u,x=±i where 



dk 



f[E^{k)]T^Jk,k + q 



{2t:Y E^Xk) - E^{k + q) + X{uj + iO) 



(8) 

Exploiting further the symmetry of the problem, in 
the limit of zero temperature we reduce the polarization 
function to the following expression 

(9) 

Here we have defined the dimensionless Fermi wave vec- 



tor vp^^ = \/ ^ — + ~ A'Cki ^^'^ the functions ^ = 
(^-fep^A ± \I^^J~-^(h^ /2aA" with 

= X cos (0k - 0q) [x cos (0k - 0q) - A^^k] ; (10) 

6,,.A = -x[(r2 + 2(A2/-a;2))cos(0k-0q) (11) 



+ sin (26*) sin (0k + 0q)] + p{Xy - 



k I 



cx = [Xy ~ x^) - x^^^ ; 
d^ = X cos (0k - 0q) - /i?k ; 



eu,A = Xy-x 



(12) 
(13) 
(14) 



jir^x 

17 



cos (0k - 0q) + sin (0k -I- 0q) sin(26')] 



Here g — m*/27r is the density of states at the Fermi 
level and we have introduced the dimensionless quan- 
tities X = q/2kF, y = oj/Asf^ V = k/kp, r — pjkp^ 
and ^k — ^{.P^^T^\i)lkp with ep = fc|,/2m* and kp = 
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FIG. 2: (a) The intra- and inter-chirality EHC in the uj — q 
plane for two different momentum orientations, 0q = 7r/4 and 
3n/4. The symbols show the plasmon dispersions (see text), 
(b) The imaginary part of the dielectric function vs energy for 
the fixed momentum magnitude q — O.lbkp and for (jiq — n/2 
and 37r/4, shown as square and round symbols, respectively. 



^/2rnFEF+p^ . The integration over v can be done an- 
alytically, yielding 
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The integration over 0^ is performed numerically. The 
derived formula (fT5|) is exact. In the limit oi 9 = (pure 
Bychkov-Rashba SOI) we recover the previous results by 
Pletyukhov and Gritsev^^, and in the limit of r = (no 
SOI) we recover the classic result by Stern^. 

For the actual calculations we use materials parameters 
suitable for InAs quantum wells with realistic values of 
the SOI parameters, r = 0.1 and 9 = 7r/8, corresponding 
to the ratio of SOI strengths a/P ~ 2.4 from^^. We take 
the electron density n = 2.55 • 10" cm^^ (Ef w 302 K) 
and the effective transverse width of the quantum well 
d = 15 nm. 

Figure [2^ shows the EHC regions and the plasmon 
dispersions for different values of the angle The 
anisotropy of the mtra-chirality EHC (the dense-hatched 
region) and of the plasmon dispersions is a small effect 
and hardly seen on the scale of figure. Meantime, the 
mter-chirality EHC is strongly anisotropic (in the long 
wavelength limit the anisotropy vanishes). Figure 
shows the imaginary part of the dielectric function vs 
energy for the fixed momentum magnitude and differ- 
ent orientations. As seen, not only the boundaries of 
EHC but also the dissipation properties within EHC are 
strongly anisotropic. In the region near the plasmon en- 
ergy, uj/AEf ~ 0.1 for q = 0.15kF, the imaginary part 
for 0q = 37r/4 is strongly suppressed with respect to its 
value for (/)q = 7r/2. 

To calculate the plasmon dispersion we solve for zeros 
of the real part of the RPA dielectric function, e{q,u!) — 
1 — v{q)ll{q,uj) where v{q) — 2tt / {nQq)F{qd) is the 
bare Coulomb interaction with kq = 14.55 the static di- 
electric constant of InAs. For the form factor F{qd) we 
use the formula (12) from Ref. [Sll, which takes into ac- 
count the transverse width d of the quantum well but 




FIG. 3: The SOI induced energy dispersion of the plasmon vs 
its propagation direction for q = O.lbkp (the left axis, square 
symbols). The parts of the spectrum which do not represent 
plasmon excitations (see text) are shown as triangle symbols. 
The dashed line plots the imaginary part of the dielectric 
function (the right axis) for q = 0.15fe_F and Aw = uj — ujq = 
1.2 K where cjq ~ 0A5Ef- 



not its asymmetric shape. The form factor has a strong 
effect on the Coulomb interaction strengtbi^ - it goes as 
1 — (1/3 — 5/4:7T'^)qd in the long wavelength limit qd — > 0, 
and as 3/(47r^g(i), in the opposite limit — > cxd. 

There are three distinct regions of EHC and the plas- 
mon disperions, as seen in Fig. [2^. In region /, which cor- 
responds to small q, the areas of inter- and intra-chirality 
subband transitions are well separated. The plasmon en- 
ergy is located within the gap between these EHC re- 
gions: these plasmons are not dampled. The plasmons 
here exhibit only a SOI induced dispersion as a function 
of its propagation orientation. At such small q, however, 
the anisotropy is not significant and eventually vanishes 
in the long wavelength limit. 

At larger values of q, in the region /// in Fig. [5^, 
the plasmon dispersion enters EHC, triggering the phe- 
nomenon known as Landau damping, i.e. decay into 
electron-hole pairs. In this regime the EHC is made up 
of several overlapping sub-regions (associated with the 
discrete quantum indices fi and A), separated by sharp 
boundaries. The imaginary part of the dielectric func- 
tion (proportional to the spectral density of electron-hole 
pairs) exhibits sharp variations across these boundaries, 
resulting in unexpectedly strong angular dependence of 
plasmon damping. In Fig. [3] we follow the evolution of 
the plasmon frequency as a function of the angle (j}q from 
to TT for q = 0.15A;f- A sharp boundary between two 
sub-regions of the EHC is crossed at (j)q ~ 7r/2. Entering 
the new region, the plasmon becomes overdamped, con- 
current with the sharp rise of Im e{q,uj), which we plot 
in the same figure on the right axis. Figure [3] shows that 
there are two ranges of directions 7r/2 ~ ~ 57r/8 and 

Ttt/S ~ ~ TT in which the plasmon cannot propagate 
due to excessive Landau damping. On the other hand, 
the plasmon is well defined around the angles = 7r/4 
and 37r/4 where the imaginary part of e(g, w) is small. 
These are the principal directions of the underlying struc- 
tural C2v symmetry. 
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FIG. 4: The dynamical structure factor vs (jiq for q — 0.15fcF. 
The sohd and dashed hues correspond, respectively, to the 
local maximum at Aui ~ 0.5 K and minimum at Alj « 1.1 K 
of the plasmon spectrum (cf. Fig. |3]) and show the single- 
peak behavior of S{q,uj). The dashed-dot-dot and dashed- 
dot lines illustrate the splitting of the structure factor peaks 
for Alj ~ 0.25 and 1.5 K. The inset shows the asymmetric 
double-peak structure of the structure factor for Acj ~ 0.7 K. 

Finally, in the intermediate // region in Fig. [2l the 
intra- and inter-chirality subbands either overlap or not 
so that the plasmon is being either damped or not de- 
pending on its propagation direction. 

In Fig. [4] we plot the dynamical structure factor 
S{q,uj) = — Im [l/e(g, w)] as a function of (fiq for uj cor- 
responding to the local maximum and minimum of the 
plasmon energy dispersion. In both cases S{q, uj) shows a 
dominant peak as a function of (f>q. As expected, the peak 
occurs at 0q = 7r/4 for uj = ujminiQ) (dashed line) and 
at (/>g = 37r/4 for lo — LUmaxil) (solid line). These peaks 
represent lightly Landau damped plasmons (the plasmon 
at (j)q = 37r/4 is less damped than the one at (j)q = 7r/4 
and therefore produces a stronger peak). The preferential 
role of these two directions comes from the C2v symmetry 
of the problem, clearly seen from the plot of the Fermi 
surface in Fig. 1. Notice that for a given uj there are 
two additional angles, at which Re s(q,uj) shows zeros. 
The structure factor, however, does not exhibit peaks at 



these angles since the large density of electron-hole pairs 
(reflected in the large value of Im e{q, uj)) overdamps the 
plasmons in these directions. These "overdamped plas- 
mons" are represented by the triangle symboles in Fig. [31 

In the range uj^axil) < uj < uJmin{<l) between the ex- 
trema of the plasmon spectrum, the height and the width 
of the peaks of S{q, uj) vs <j)q show a smooth evolution: 
with increasing uj one peak diminishes, the other grows, 
and vice versa. Thus, in this intermediate region the 
structure factor has two peaks, located at 0q = 7r/4 and 
37r/4, which constitute an asymmetric doublet, shown in 
the inset of Fig. [4] In the energy regions above the mini- 
mum or below the maximum of the plasmon spectrum at 
given q (i.e., for uj > ujminiq) or w < ujmaxiq)), Re e{q,uj) 
vs (jjq shows two zeros around 7r/4 or 37r/4 so that each 
peak of S{q, uj) splits into two peaks located symmetri- 
cally above and below the angle (jjq — Tr/A (the dash-dot 
line) or 37r/4 (the dash-dot-dot line). 

In the case of a = ±/3 (see Fig. [T|d) the linear spin- 
orbit couplings do not affect the plasmon spectrum: the 
plasmon damping vanishes and the structure factor is a 
delta-function for all momentum orientations. For this 
special case there is a global spin quantization axis - one 
of the principal C2v axes - so that the electron gas is 
split into two uncoupled spin components, whose circu- 
lar Fermi contours are shifted from the origin in oppo- 
site directions. Each component gives an isotropic col- 
lective response, as guaranteed by Galilean invariance. 
Cubic spin-orbit terms, which spoil this effect, are typi- 
cally much weaker in quantum wells. 

In conclusion, we have shown that plasmon dynamics 
(spectrum and damping) is strongly anisotropic in real- 
istic zinc-blende quantum wells, due to the interplay of 
two different forms of spin-orbit interaction. Experimen- 
tal observation of this anisotropy would be of fundamen- 
tal interest and could open the way to new techniques for 
controlled directional plasmon filtering potentially useful 
both for spintronic and plasmonic devices. 
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